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Abstract. We study in detail the generation and relaxation of quantum coherences 
(entanglement) in a system of coupled polariton traps. By exploiting a Lie algebraic based 
super-operator technique we provide an analytical exact solution for the Markovian dissipative 
dynamics (Master equation) of such system which is valid for arbitrary cluster size, polariton- 
polariton interaction strength, temperature and initial state. Based on the exact solution of 
the Master equation a,t T — OK, we discuss how dissipation affects the quantum entanglement 
dynamics of coupled polariton systems. 



1. Introduction 

Coupled bosonic systems appear naturally in a great variety of physical situations ranging 
from quantum optics, atomic/molecular physics to condensed matter physics. In the last field 
exciton/polariton systems are among the most notable examples. Similarly to ultracold boson 
atoms confined in optical lattices, trapped exciton/polariton systems are predicted to present 
different collective phases such as Mott insulator or superfluid phases. Experimental control 
of polariton trapping is rapidly improving. Among the most promising candidates are stress 
induced trap systems [1] and metal thin-film deposition on microcavities [2]. A proposal for the 
generation of polarization entangled photon-pairs from electric gate controlled polariton trapped 
systems has been recently put forward \3\- However, for full exploitation of the capabilities of 
semiconductor microcavities based entanglement, the polariton entangling power itself requires 
complete characterization. In the present work, we focus on polariton trapped clusters as open 
quantum systems evolving under nonequilibrium or dissipative conditions. The time evolution 
of the entanglement generated by polariton-polariton interactions jointly with dissipative effects 
is considered. 

In order to analyze the dissipative quantum dynamics of a trapped polariton system we use 
a Born-Markov reduced dynamics approach. The solution of the resulting Master (Lindblad) 
equation (ME) describes the time evolution of the polariton system of interest. While deviations 
from a strictly Markovian behavior may occur in real microcavities [1], our analytical results 
may be used as a first step in exploring more complex relaxation dynamics in multitrap polariton 
systems. We restrict ourselves to the case where the system is deep in the Mott insulator or fully 
localized regime where polariton hopping is neglected but polariton-polariton interactions remain 
important. To find the full density matrix of the coupled polariton system, we use an exact Lie 



algebraic superoperator method for solving the ME, generalizing previous similar treatments 
[U Ej to arbitrary number of coupled systems. By applying these results to a double trap, we 
focus our attention to noise related properties of the entanglement creation and preservation 
in such condensed matter highly nonlinear systems. We quantify the amount of entanglement 
in terms of the logarithmic negativity which has been proven to be valid as an entanglement 
measure for any mixed state of arbitrary bipartite systems [7]. Furthermore, this entanglement 
measure is easily computable having additionally an operational interpretation. We found that in 
the Mott regime, the polariton-polariton interaction generates naturally entanglement which is 
degraded by dissipation on a shorter time scale as compared with the population decay, and with 
an entanglement maximum being an increasing function of the initial number of polaritons. The 
entanglement duration time shortens as the polariton decay rates or the initial trap populations 
become large. 



2. Theoretical method 

Within the Born-Markov approach, the ME for the reduced system's density operator p can be 
written as (?i = 1) 
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dt 
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with the first term on the right hand side describing the coherent or unitary evolution term, 
whereas the second term, or Lindblad term, is associated to the noisy dynamics induced by the 
reservoirs. Thus, we address the study of the relaxation dynamics of M coupled polariton traps 
described by the Hamiltonian 
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where (a„) are boson creation (annihilation) operators in trap n. The coupling matrix 
A = {A„.m} is a real symmetric matrix (hereafter diagonal terms will be denoted simply as 
An,n = Lindblad terms, as taken to describe reservoirs acting independently on each 

polariton trap, are given by 
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where coefficients -D„ and Pn are related to the dissipation and pumping mechanisms, 
respectively, for the n-th trap. In this equation, there are no cross terms of the type ci^(iyYi 
because we consider that decay and pumping on one trap is independent from other traps. 

In order to evaluate the time evolution of the polariton traps density operator we proceed as 
follows. We define super-operators in the following form 

Ki")o = aidan , K^^'^d = andai , ^'^"^0 = -^(4o„6 + O^On + O) (4) 

and AfnO = [aJ^On, O], where O denotes any polariton operator. Furthermore, given the fact that 
[djjd„d]„am, /5(t)] = — Mn (-^o™'' + ^) + J^m (^Kq^^ + p{t), the ME for the coupled polariton 
system, Eq.([l|, can be written as 
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with U^n = LOn " A„ - Em^^n ^n.m and (in = Dn + Pn + i^nMn + iY.m^n^n,mJ^m- As the 

super-operators Mn commute with any other super-operator, i.e. [MrmJ^n] = [■f^m,K^'^] = 
[Mm, k'q''^] = 0, the solution to Eq.(l5|) has the form 
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where p{0) denotes the initial state. It is straightforward to demonstrate that the superoperators 
defined by Eq.(j3]) close in the su(l,l) Lie algebra 
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By using the su(l,l) Lie algebra structure of the superoperators as given in Eq.([7]), the ME 
solution expressed by Eq.Q, can finally be written in the full disentangled form 
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where the temporal coefficients c^^\t) and Cg"'''(i) are functions of Q.n (details are presented 
elsewhere In this way any expectation value of polariton observables could be analytically 
evaluated. 



3. Entanglement nonequilibrium dynamics 

Given the general solution of the ME in Eq.([8]) an analytical expression can be obtained for the 
cluster density operator [8]. However, in the present work we focus on the zero temperature 
case for a two trap polariton cluster (M = 2). The same algebraic results could also be used 
to model a single trap with two polariton species, as it could be the case for a system formed 
by polaritons of different spin polarizations confined in the same spatial region. We consider 
that the system is initially prepared, by a resonant short laser pulse, in a purely separate 
coherent state, i.e.. /5(0) = |ai,a2 >< ai,a2|, where \aj > denotes a coherent state in the j-th 
trap, j = 1,2. The initial mean number of polaritons in trap j is given by |ajp. After this 
initialization step, interacting polaritons evolve in time in presence of dissipation while, in the 
vanishing temperature limit, incoherent pumping effects are not present {Pn = 0). Moreover, 
we take identical dissipation decay rates Dn = T. We take as the unit of energy the inter-trap 
interaction strength Ac = 1 and consequently time is measured in units of A~^. 

Now the question is: how to identify entanglement between two bosonic systems?. By 
contrast with two qubit systems for which entanglement can be unambiguously identified by 
the concurrence |9j, for two continuous variables as it is the case here, sufficient and necessary 
conditions for detecting entanglement have only been established for the special case of Gaussian 
states \10\ lllj. However, in general the polariton trap state is a mixed non-Gaussian state. 

We quantify the polariton entanglement by computing the logarithmic negativity LN(p) = 
log2\\p^'^\\, where is the trace norm of the partial transpose density operator [7J. This 

latter quantity can be easily evaluated as ||/0''^'^|| = 2N{p) + 1 where N{p), or negativity, is the 
absolute value of the sum of the negative eigenvalues of p'^"^. Based on the exact analytical 
solution for the ME given in Eq.Q, the required transpose operation on the system's density 
operator is trivial for any time from which the negative eigenvalues of are numerically 
evaluated. For any separable or unentangled state LN{p) = N{p) = 0. Obviously, the 
entanglement fully vanishes for Ac = 0. To focus the discussion and for demonstrational purposes 
we consider here the case of identical traps, i.e. uji = L02 and Ai = A2, which are also identically 
prepared ai = 02- 
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Figure 1. Logarithmic negativity 
as a function of time for different 
polariton self-interaction strengths 
and dissipation rates. See text for 
parameters. 

In Figure (1) we plot the time evolution of the polariton entanglement, LN[p{t)], for a two 
trap cluster initially prepared with |aj|^ = 6 (the amount of entanglement does not depend on 
the frequency ujj of each trap). Four different values of dissipation rates F = 0.01,0.1,1 and 
3 are considered and correspond in Fig.(l) to colors red, blue, black and green, respectively. 
Additionally, for each decay rate, results for two extreme values of the ratio between the intra- 
trap and inter-trap interactions are displayed: Xj/Xc = 1 solid lines and Xj/Xc = 10 dashed lines 
(colors are still as described above). The first observation is that entanglement is immediately 
generated by the polariton-polariton interaction, Ac, and it is preserved for a significant fraction 
of the time that polariton population in the traps is maintained. The cluster system starts its 
evolution in a separable state rising its entanglement degree in an almost linear way with a slope 
depending directly on the initial polariton population, a clear indication of how important is the 
effective non-linear polariton interactions for generating entanglement. At long times the system 
ends again in a separable state, as it should be because the traps become empty and thus the final 
state is nothing but the double vacuum state. At very low decay rates entanglement oscillations 
are visible due to the quasi-unitary evolution of inter-trap correlations governed exclusively by 
Ac [H] • Note that in this limit of weak dissipation the entanglement dynamics does not depend on 
the intra-trap polariton interactions Xj, as it is evident from the almost identical solid and dashed 
red curves in Fig.(l). As the dissipation becomes important populations and the entanglement 
decay faster. The polariton population follows a very simple exponential decay form given by 
< OjCij > (t) = |ajpe~^'"* while for the entanglement a more pronounced dependence on the 
polariton self-interactions are evident. Thus, dissipation and polariton self-interactions jointly 
degrade significantly both the maximum attainable entanglement and the time interval where 
entanglement is visible. 

Figure (2) shows the time evolution of LN(p) for different polariton initial populations 
|ajp = 10 (red curve), |ajp = 8 (blue curve), |ajP = 6 (black curve) and |ajp = 3 (green 
curve). Other parameters are: Ai = A2 = — 1 and F = 1. It is evident that the rising time of the 
entanglement is shorter as the polariton population increases, emphasizing previous comments on 
the importance of inter-trap polariton-polariton interactions for an efficient and rapid generation 
of entanglement. However, on the other hand, the temporal width of the entanglement curve 
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Figure 2. Logarithmic negativity 
as a function of time for different 
polariton populations. See text for 
parameters. 



decreases as a function of the polariton population, indicating the high level of non-linear effects 
occurring in the present system. 

Some numerical estimates are in order. According to Ref.[2] typical engineered trap radius 
range from 5 to 50 /im. Pumping pulse lasers of low intensity on the order of 1 mW give rise to 
a polariton density of 10* — 10^*^ cm^^. Thus, for experimentally achievable trap sizes, typically 
10 — 100 polaritons could be found in a single small trap. Our unit of time is fixed by the 
inter-trap (or inter-species) interaction strength Ac- According to Ref.[3j, Ac ~ 10 — 100 fieV, 
which for Xct ~ 1 yields to f ~ 40 psec. Hence, the theoretical results we discuss above are in 
the lower range of experimentally accessible polariton densities and time resolution, requiring 
very clean microcavities samples and high quality factors Q, for entanglement detection to be 
observable. However, our analytical results should be indeed valuable as a starting point for 
more numerically involved calculations in realistic situations. 

4. Conclusions 

By solving exactly the Master equation for an interacting trapped polariton system we have 
found that polariton entanglement can be immediately produced from an initial separable 
multi-coherent state. The entanglement starts rising linearly with a slope depending on the 
initial trap populations. The entanglement maximum is an increasing function of the initial 
number of polaritons. However, the entanglement duration time shortens as the initial trap 
population becomes large. In any case polariton entanglement observation would require very 
good microcavities with weak dissipation rates. 
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